Dynamics in inhomogeneous liquids and glasses via the test particle limit 
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We show that one may view the self and the distinct part of the van Hove dynamic correlation 
function of a simple fluid as the one-body density distributions of a binary mixture that evolve in time 
according to dynamical density functional theory. For a test case of soft core Brownian particles the 
theory yields results for the van Hove function that agree quantitatively with those of our Brownian 
dynamics computer simulations. At sufficiently high densities the free energy landscape underlying 
the dynamics exhibits a barrier as a function of the mean particle displacement, shedding new 
light on the nature of glass formation. For hard spheres confined between parallel planar walls 
the barrier height oscillates in-phase with the local density, implying that the mobility is maximal 
between layers, which should be experimentally observable in confined colloidal dispersions. 

PACS numbers: 61.20.Lc, 82.70.Dd, 64.70.Pf 



The van Hove function G(r, t) for the probability of 
finding a particle at time t at a distance r from the ori- 
gin, given that there was a particle at the origin at time 
t = 0, characterizes dynamical phenomena in Condensed 
Matter on a nanoscopic scale Q, Q- Recent measure- 
ments of G(r, t) using confocal microscopy were aimed at 
investigating e.g. dynamical heterogeneities in colloidal 
hard sphere suspensions 3] , or the devitrification of col- 
loidal glasses upon adding nonadsorbing polymers 
The Fourier transform of G(r,t), the intermediate scat- 
tering function, can be measured with inelastic scattering 
techniques |2|. 

There is no consensus on a theoretical approach for 
investigating the microscopic dynamics in dense inhomo- 
geneous liquids, relevant to studying e.g. the glass tran- 
sition for liquids adsorbed in nanoporous materials 0. 
This situation contrasts with the static case, where mod- 
ern density functional theory (DFT) provides a powerful 
means for investigating the significant effects on phase 
behaviour and structural correlations that spatial con- 
finement and external fields may induce 0, El- DFT 
operates on the level of the one-body density distribu- 
tion p(r), where r is the space coordinate. For a range 
of situations the time evolution of the density profile, 
p(r, t) , as induced by a time-dependent external poten- 
tial, V ex t(r,t), has been shown to be well described by 
dynamical density functional theory (DDFT) 0,0, see 
e.g. the successful applications in Refs. 0>ilE3- How- 
ever, currently no similar theoretical framework exists for 
calculating dynamical two-body correlation functions, as 
required to address problems such as those mentioned 
above. 

In the static case a close relationship between p(r) and 
the pair distribution function, g(r), describing the prob- 
ability of finding a pair of particles separated by a dis- 
tance r, is established by Percus' famous test particle 
limit pg(r) is the one-body density distribution of 



a fluid exposed to the influence of an external potential 
V cx t( r ), that represents a 'test particle' fixed at the origin, 
given by V ex t( r ) = V( r )i where V(r) is the interparticle 
pair potential, and p is the bulk density. 

In this Letter we generalize the test particle limit to 
dynamical correlation functions allowing us to use DDFT 
to calculate the van Hove function in bulk and in inho- 
mogeneous systems. We test this approach by comparing 
results to those from our Brownian dynamics computer 
simulations for a simple Gaussian core model (GCM) for 
a macromolecular solution in bulk, and indeed find very 
good agreement. For hard spheres we estimate the glass 
transition to be at bulk densities pa 3 ~ 0.9, where a 
is the hard sphere diameter. Upon confining the hard 
sphere system between parallel hard walls, we find the 
local mobility to be maximal where the local density is 
minimal. Our approach allows for inspection of a well- 
defined underlying free energy landscape, which requires 
solving a corresponding equilibrium DFT rather than the 
foil DDFT, easing numerical burdens. 

Consider a fluid of N particles that interact via a pair 
potential V(r). The system may be viewed as a binary 
mixture of species s and d, in which two of the pair po- 
tentials Vij(r) between particles of species i,j = s,d are 
equal, V dd (r) = V sd (r) = V(r), and V ss (r) = 0. Fur- 
thermore species s consists of only a single (test) par- 
ticle, located at position ro at t = 0, while species d 
consists of the remaining N — 1 particles of the system. 
Our aim is to relate the one-body density distributions 
of this mixture, p s (r,t) and p d (r,t), to the self and dis- 
tinct part of the van Hove function of the pure fluid, 
G„(r ,r,t), and G d (r ,r, t), respectively. [Recall that 
G(r,t) — G s (r,t) + G d (r,t).] We must choose suitable 
initial conditions for p s (r,t — 0) and pd(^,t = 0). In the 
case of a bulk system let the one-body density distribu- 
tion of species s at time t = be p s (r,t = 0) = S(r), 
where S(-) is the Dirac delta function and we have cho- 
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sen the coordinate system such that ro = 0. Species 
d is initially at equilibrium, i.e. pd(r, t = 0) = pg{r). 
For subsequent times, t > 0, we make the identification 
G s (r,t) = p s (r, i) and G d (r,t) = p d {r,t). 

There are two ways to proceed. One is to carry out 
computer simulations on the level of particle coordinates 
(detailed below). The alternative is to use a theory that 
operates on the one-body level, and DDFT is an obvious 
choice. Generally there are only formal results for the ex- 
act equations for the time evolution and hence one must 
resort to approximations 0J1,H>E1> sucn as the theory 
of Marconi and Tarazona [7j for Brownian dynamics: 

dp, ( 5F[p s ,p d }\ . 

where Y is a friction coefficient characterising the drag of 
the solvent on the particles, and F[p s , p d ] is taken to be 
the equilibrium Helmholtz free energy functional: 

F[ Ps , Pd ] =k B T J drp«(r)[ln (p 4 (r)A 3 ) - 1] (2) 

2=5, d 

+ F cx [p s ,p d ] + J drV r oxt (r)[ / o s (r) + p d (r)], 

where the first term on the r.h.s. is the Helmholtz free 
energy functional of the (binary) ideal gas, A is the (ir- 
relevant) thermal wavelength, T is the temperature, k B 
is the Boltzmann constant, and F ex [p S) pd] is the excess 
contribution to the Helmholtz free energy due to inter- 
actions between the particles [l4|. 

As a test case we consider the GCM, that describes 
the soft interactions of polymer coils, star polymers, or 
dendrimers in solution 15], and is defined by V(r) = 
eexp(— r 2 /i? 2 ); R is the radius of the particles and e is 
the energy cost for complete overlap of a pair of particles. 
We use a simple mean field approximation, F ex [p s ,p d ] = 
\ J2i,j= s ,d I drdr 1 p % {r)Vij{\v - r'\)pJr'), which has been 
shown to be reliable for this model |9( , and integrate Eq. 
l|T)l numerically with the given initial conditions forward 
in time. In Fig. ^ results are displayed for both parts of 
the van Hove function at four different times for a typ- 
ical statepoint. We also display results from Brownian 
dynamics (BD) simulations of 300 particles, with a time 
step of 10 _3 r B , where t b = R 2 /(Yk B T). The very good 
agreement between the DDFT and the BD simulation re- 
sults attests to the validity of our scheme. Nevertheless, a 
remark about the more general applicability of the DDFT 
in Eq. Q is in order. The relevant dynamical variable is 
the root mean square particle displacement w(t), defined 
as the width of G a (r, t) [and hence of p a (r, t)} via 

poo 

[ w (t)} 2 =4tt / drr 4 G s {r,t). (3) 
Jo 

For t — > oo, w(t) 2 = 6Drt, where Dl is the long time 
self diffusion coefficient 2}. However, for not too low 




FIG. 1: The self part G s (r,t) (dashed lines, squares) and 
scaled distinct part Gd(r,t)/p (full lines, circles) of the van 
Hove correlation function for the GCM fluid with bulk den- 
sity pR 3 = 0.2263 and k B T/e = 0.5, at times t* = tf tb = 
0.01,0.1,0.2 and 0.4 (a,b,c,d, respectively), as a function of 
the (scaled) distance r/R. Shown are results from DDFT 
(lines) and BD simulations (symbols). 

temperatures and t — » oo, Eq. predicts that w(t) 2 = 
6Dst, where D$ = k B T/Y is the short time self-diffusion 
coefficient, and for the GCM Dj, ~ D$ 13]. However, 
for particles with a hard core Dj, < D$ and the DDFT 
in Eq. Q is not reliable for t ^S> t b . To circumvent this 
problem one could use a DDFT such as that of Ref. [13j , 
which guarantees the correct long time behaviour - see 
also the approach of Ref. ^(| . 

Here we choose a different direction and investi- 
gate the free energy landscape underlying Eq. 
in more detail. We consider the 'fruit-fly' of liq- 
uid state physics, namely a one-component sys- 
tem of hard spheres. This is described by the 
Ramakrishnan-Yussouff (RY) approximation ll^j : 
F cx [ Ps ,p d ] = Vf cx ( P ) + Ei= s ,dfL(p)JdrApi(T) - 
lE ili=s ,J drdr'Ap, (r)cy (|r - r' | )Ap j (r') , where / cx (p) 
is the bulk excess free energy per volume V, f' cx {p) — 
df ex {p)/dp 1 Ap t (r) = pi(r)-pi, with p d = p, p s = 0, and 
the Percus-Yevick approximation [2j for the pair direct 
correlation function c i3 (r) ; we set c ss (r) = to model the 
absence of test particle self-interactions. Chosen here for 
its simplicity, this approximation for F ex [p a ,p d ] is crude, 
but sufficient to explore qualitative behaviour. We use 
equilibrium DFT to elucidate a pathway that closely re- 
sembles that of the full DDFT by minimising the free en- 
ergy functional |(J, i.e. solve SF/Spi(r) = pt for i = s,d 
where pi are Lagrange multipliers to fulfill J drp s (r) = 1 
and / dr(pd(r) — p) = p J dr(g(r) — 1) = const. As a 
constraint we control w via an associated Lagrange multi- 
plier, A, which is formally equivalent to treating an auxil- 
iary external field acting on component s with the simple 
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harmonic form V}^(r) = Xr 2 /a 2 . 

By calculating F as a function of A (and hence w) 
for a given statepoint, we may determine the free en- 
ergy landscape that governs the time evolution of G(r, t). 
We find that at low densities F(w) decreases monoton- 
ically with w, but as p is increased above a threshold, 
p > p* , F(w) becomes non-monotonic and develops a 
barrier; the barrier height AF, i.e. the difference between 
the minimum and maximum, grows upon increasing p, 
see Fig. [JJ While for both small and large values of 
w, p s (r;X) is nearly Gaussian, for intermediate values, 
particularly when p > p* , p s (r; A) has pronounced non- 
Gaussian characteristics; Fig. [2K shows results for p s (r) 
for three typical values of w. For states that correspond 
to barrier crossing p s (r) exhibits a pronounced shoulder, 
taken as a signature of dynamical heterogeneity 0, 0] . 
However, assuming the Gaussian form for all values of 
w, as in the main plot of Fig. [JJ does not affect F(w) 
strongly, and only slightly overestimates AF (see Fig. [SJd 
for a comparison of results from free minimisation and 
from the Gaussian parametrisation). The emergence of 
a barrier in F(w) leads to a trapping in the DDFT for 
p > p* whereby, as t — > oo, w(t) — > wq < oo, suggest- 
ing that the system is non-ergodic. However, in reality, 
as long as the barrier is sufficiently small, it can eventu- 
ally be overcome. Our expression for F cx is approximate 
and neglects some contributions to the free energy due 
to fluctuations 0, 0] , so the particular DDFT in Eq. (JIJ 
does not describe the barrier crossing; one could include 
a stochastic noise term, following Ref. 0], to take this 
into account. 

We may use our results to estimate when the hard 
sphere fluid becomes trapped in a glassy state by con- 
sidering the time t it takes to cross the barrier, which 
scales as r - e AF/fcBT [H H|. We estimate the lo- 
cation of the glass transition by determining the state 
point at which AF reaches a particular value; as a rough 
criterion we take lOfcsT. From the results in Fig. [21 one 
would expect this to occur around p g <r 3 — 0.9, which 
is somewhat below the result p g cr 3 — 1.1 found exper- 
imentally for colloidal hard spheres j2l|. We attribute 
this under-estimation to our choice of simple (RY) free 
energy functional; use of a more reliable functional [22 ] 
may improve results. 

We emphasise that our approach is distinct from theo- 
ries describing the glass transition as freezing into an ape- 
riodic lattice |2jj ; we do not find a thermodynamic phase 
transition. There is some similarity between our phi- 
losophy and that underlying the successful 'non-classical 
theory of nucleation', where DFT is used to calculate the 
free energy barrier to nucleating a liquid droplet in the 
oversaturated gas |24| . The results presented in Fig. [2J 
bear much similarity to those pres ented in recent work 
by Schweizer and Saltzman (lfl l20j| who also construct a 
theory for the free energy of a single particle in the cage 
of its neighbours. One key difference is that we calculate 




FIG. 2: The scaled free energy F/ksT as a function of the 
root mean square particle displacement w [see Eq. 0], cal- 
culated with the RY functional for hard spheres at densities 
per 3 = 0.6,0.7,0.8 and 0.9 (as indicated by the arrow) and 
shifted such that F(w = 3) = 0. (a) Self part of the van 
Hove function on a logarithmic scale as a function of r/a for 
pa 3 = 0.78 and three typical widths w. Note the pronounced 
non-Gaussian character for w/a = 0.45. (b) Comparison of 
results for F/ksT as a function of w/a from free minimisa- 
tion (dot-dashed line) and a Gaussian parametrisation for G s 
(solid line) for per 3 = 0.78. Also shown is the free energy bar- 
rier AF between minimum and maximum (arrow) and the 
values of w/a (circles) for which results are shown in (a). 



the free energy of the entire system, including contribu- 
tions from Pd{r) and p s (^), whereas their theory is for 
the free energy of the test particle only. Furthermore our 
approach is readily applicable to inhomogeneous systems. 

We have investigated how the dynamics of hard spheres 
is affected by confinement between two planar parallel 
hard walls separated by a distance H, described by an 
external potential V ox t(z) that vanishes for —H/2 < z < 
H/2, and is infinite otherwise; z is the coordinate per- 
pendicular to the walls. The presence of the walls causes 
the one-body density, p(z), to be oscillatory with maxi- 
mal value at contact with the walls; see Fig.[3Ji for results 
for a system with H = 4er in chemical equilibrium with 
a bulk of density per 3 = 0.8. Due to symmetry the van 
Hove function depends on the initial position zq of the 
test particle, but not on its initial lateral coordinates, 
taken as Xo = yo — 0. We have used equilibrium DFT to 
determine Pd(r), assuming p s (r) to be a Gaussian (recall 
that in bulk this assumption made little difference to the 
value of F). As an example we show in Fig. [31) results for 
p s (r) and Pd(r) for the case where the test particle was 
initially in the midplane, zq = 0. We choose conditions 
(see below) where p s (r) is still peaked around r = 0; 
Pd exhibits structure due to packing effects caused by 
the walls, as well as by packing around the test particle. 
Quite striking is the appearance of the hexagonal shape 
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FIG. 3: Correlation functions for hard spheres confined be- 
tween parallel plates of separation distance H — 4a. (a) Den- 
sity profile p(z) (full line) as a function of z/a and free energy 
barrier AF(zq) (symbols; the dotted line is to guide the eye). 
Also shown is the corresponding bulk value of AF (dashed 
line), (b) The (scaled) self part p 3 (r)a 3 /25 and distinct part 
p,j(r) of the van Hove function as a function of the scaled 
coordinates perpendicular and parallel to the wall, z/a and 
x/a, respectively, evaluated at the minimum of the free energy 



landscape for bulk density pa 



0.8. 



of the first coordination shell; for large planar distances 
Pd(z, x — > oo, y = 0) = p{z). Similar to the bulk case we 
find a free energy barrier AF that now depends on zq (the 
profiles in Fig. 03 are shown at the minimum of F.) Fig. 

shows strong oscillations of AF(zq) around its value 
in bulk; the oscillations are in phase with the oscillations 
of p(z). Since the barrier hopping time r ~ e AF/k B T ^ 
the particles in the layers (i.e. where p(zo) is large) are 
less mobile than those between the layers. For the case 
in Fig. [3t> the particles at z = ±0.5cx have a mobility 
(~ I/r) that about 50% of the bulk value. This predic- 
tion can be tested by experiment with confined colloidal 
hard spheres or with computer simulations. 

In conclusion we have presented a general framework 
for calculating the van Hove two body correlation func- 
tion both in bulk and in inhomogeneous fluids. We have 
proposed a method for determining the free energy land- 
scape for a particle in the cage of its neighbours and 



have demonstrated the applicability to inhomogeneous 
systems, where e.g. mode coupling theory cannot be 
easily applied |l3| . 
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